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ABSTRACT 

We carry out a Bayesian model selection analysis of different dark energy 
parametrizations using the recent luminosity distance data of high redshift su- 
pernovae from Riess et al. 2007 and from the new ESSENCE Supernova Survey. 
Including complementary cosmological datasets, we found substantial evidence 
(Aln(E') ~ 1) against a time-varying dark energy equation of state parameter, 
and against phantom dark energy models. We find a small preference for a stan- 
dard cosmological constant over accelerating non-phantom models where w is 
constant, but allowed to vary in the range —1 to —0.33. 

Key words: cosmology, dark energy. 



1 INTRODUCTION 



Over the last few years, observations of luminosity dis- 
tances of Type la supernovae (SN-Ia) have established 
that the expansion of the Universe is a cceler a ting (see e.g 



that the expansion oi the Universe is accelerating (s ee e.g. 
IRiess et al . (1998), IPerlmutter et al.1 (Il999h. IRiess et al.1 



(|2004l ). lAstier et al.1 (|2005h . iBassett et all (|2004h ). This 
result is now well confirmed and complemented by a 
large amount of independent observations such as, for 
example, the angular- diameter distance vs. redshift rela- 
tion measure d by Baryonic A c oustic Oscillations (BAO) 
experiments lEisenstein et al.l ([2005) , the distortion of 
backg round images mea sured by weak lensing experi- 
ments |jarvis^£^l] (|20Q5h , the distance to the last scatter- 
ing surface measured by Cosmic M i crowa ve Background 
(CMB) experiments ISpergel et al ] (|2006l), galaxy clus- 
tering (Large Scale St r uctur e, LSS) (see lEfstathiou et al.l 
(|2002h : lTegmark et al.1 (120061 )) and, finally, the Integrated 
Sachs Wolfe effect, correl ating LSS with CMB (see e.g 
Giannan tonio et al 1 (l2006h ). 

The recent analysis of Rie ss et all (|2007l ) has fur- 
ther confirmed in an impressive way these results, re- 
porting the discovery of 21 new SN-Ia with the Hub- 
ble Space Telescope (HST). Together with a recalibra- 
tion of previous HST-discovered SN-Ia, the full sample of 
23 SN-Ia at z > 1 provides the highest-redshift sample 
known. This dataset has then been analyzed in combina- 
tion with some of the aforementioned datasets providing 
new constraints on severa l dark e nergy properties (see e.g . 
Riess et al.1 (1200 /t) . lAlaml (120061) , [Gong and Wan j tOOdi ) 
Nesseris &; Per i volar op oulosl (|2006r ) jBarger et al.l (|2006r )). 

This increasing quality and number of experimen- 



tal datasets is finally opening the possibility of falsify- 
ing cosmological theories and of discriminating between 
different theories. There have been many proposed ex- 
planations for this acceleration: the Einstein's cosmolog- 
ical constant, a new fluid with negati ve pressure (con- 
stant or varying with time) (see e.g. IPeebles &; Ratral 
(2003)) or a modification of general rela tivity (see for 
example ID vali, Gabadadze, Porratl (|20QQh ) . However, to 
date, none of them is supported by a well-established fun- 
damental theory. Moreover, since we know (almost) noth- 
ing about dark energy there is in principle no theoretical 
limit to the number of parameters that one might use to 
characterize it. 

It is therefore timely not only to constrain the param- 
eters of a specific dark energy model but also to establish 
reliable criteria to choose between different models. 

As po inted out in iMukheriee, Parkinson, Li ddle 

(2006) and I Liddle et al ] (|2006aT ). there is an important 
difference between parameter fitting and model selection. 
In the first case we work in the context of a single theo- 
retical model framework, and establish which choice of 
model parameters gives the best fit to the data. In a 
Bayesian interpretation, with flat priors the most likely 
model parameters will simply have the maximum likeli- 
hood and the lowest x 2 regardless of the degrees of free- 
dom. More complicated models (with a larger number of 
free parameters) will normally produce better fits. How- 
ever, in model selection, we wish to know which model 
is favoured, regardless of the values of the parameters. 
In choosing models, we should also look for simplicity, 
following somewhat a principle based on Occam's Razor. 
Simplicity is of course in contrast with more parameters 
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and we need then to weight the "need" for extra param- 
eters. This can be accomplished by using the Bayesian 
Evidence, defined as the probability of the model given 
the data, and given by the average likelihood of a model 
over its prior parameter space: 



E = P(D\H) 



-I 



d0P(D\e,H)P(9,H) 



(1) 



where H is the model considered, is the vector of the 
model parameters, D is the data, P(0, H) is the prior and 
P(D\0, H) is the likelihood. For flat priors on the models, 
the probability of the model given the data is proportional 
to the Evidence. 

In general, comparing two models, the term 
P(D\0,H) will be larger for the more complicated one 
(which has more parameters or the same number of pa- 
rameters but with larger prior space), but, at the same 
time, it will have a lower P(9, H)dO compared to the sim- 
pler model. The first term indicates how well the model 
fits data, the second one indicates how simple is the 
model. Following the equation above we can assume that 
the best model will be the one with the greatest Evi- 
den ce 

|jeffrevsl(|l96lh provides a useful guide to discriminate 
the difference between two models with E\ and E2: 



1 < A \n(E) < 2.5 (substantial) 
2.5 < A \n(E) < 5 (strong) 

5 < A \n(E) (decisive). 



(2) 
(3) 

(4) 



The Bayesian Evidence can therefore help us to choose 
between models because it establish a tension between 
the simplicity of a model and its power of fitting data. 

In this paper we aim exactly to make use of cos- 
mological model selection methods in order to discrim- 
inate between dark energy models. In the next sections 
we therefore compute Bayesian Evidence for a large set 
of models and we compare this value with the one ob- 
tained for the concordance model: a flat universe with 
cold matter and a cosmological constant. More specifi- 
cally, in Section 2 we present our analysis method and 
the dark energy models considered. In Section 3 we show 
the results of our analysis and we derive our conclusions 
in Section 4. Our paper follows the research lines already 



i nvestigated by previous papers (see e.g. iLiddle et al 



(1200601). ISaini et al.1 ri2Q04MElgarov fe MultamakH <|2006h 
and ISzvdlowski et al.l ([20061 ) ) which we will complement 
and extend by using a new and independent algorithm 
for computing evidence (which we will illustrate in Sec- 
tion 2.3), a larger set of dark energy models and finally 
by considering more recent SN-Ia datasets. 



1 We want to stress that the "best" model might not be in 
general the "true model" ; Nature decides what is true, not our 
personal aesthetic behaviour in simplicity. 



2 ANALYSIS METHOD 



2.1 Theoretical Framework 



We restrict our analysis to flat, Friedmann-Lemaitre uni- 
verses, with the redshift evolution of the expansion rate 
given by: 



H 2 (z) = Hi 



n m (i + zf + (1 - n m ) 



Px(z) 



MO) 



(5) 



where Q m is the energy density parameter in matter, H(z) 
is the Hubble constant and px is the dark energy density 
given by: 



Px(z) 
Px(0) 



(6) 



where w(z) is the equation of state parameter defined 
as the ratio of pressure over density of the dark energy 
component Px — w(z)pxc 2 . 

We consider different parametrizations of the dark 
energy equation of state parameter w(z). The simplest 
model is the usual flat ACDM model with fixed equation 
of state w = — 1 (MODEL I). We then let w vary, assum- 
ing it is small enough to lead to acceleration. MODEL 
II has constant w with a flat prior in the range — 1 ^ 
w ^ —0.33. In MODEL III, we expand the prior range 
to allow phantom dark energy models, constant w with a 
flat prior — 2 ^ w ^ —0.33. We also consider dynamical 
dark energy models where w can depend on redshift. In 
particular we co nsider a linear dependence on scale factor 
a = (1 + z)- 1 as lChevallier fe Polarskil (|200lh : 



w(a) = Wq + w a (l 



(7) 



with -2 < w ^ -0.33 and -1.33 ^ w a ^ 1.33 
(MODEL IV). The above model is a low redshift ap- 
proximation that may break at higher redshift. In this 
respect it is useful to include a more sophisticated 
parametrization that takes in to account the high red- 
shift behaviour. We consider two possibilities. The first 
one is the one proposed by Hannestad and Mortsell (see 
lHannestad fc Mortselll (|2004)), where: 



w(a) — wqWi 



a q + a q s 



w\a q + woa q 



(8) 



In this model (MODEL V), the equation of state changes 
from wo to w\ around redshift z s — 1 — l/a s with a 
gradient transition given by q. The priors are flat within 
-2 < wo < 0, -2 ^ wi < and < q ^ 10. 

The second one i s the parametrization introduced by 
lUpadhve et~ai1 (|2005h where 



w(z) = Wo + W\Z 
for z < 1 and 
w(z) = Wq + Wi 



(9) 



(10) 



for z ^ 1 (MODEL VI). In this case we choose flat priors 
in —2 ^ wo ^ —0.2 and —4 ^ w\ ^ 2. In analyzing each 
model, the priors for the other parameters are flat within 
the ranges 
0.1 ^ On ^ 0.5 

56^ HoKkms- 1 Mpc' 1 ) ^ 72. 
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2.2 Cosmological Datasets 

The dark energy models are then com pared with the data 
follow in g the approa c h desc ribed in IWang &; Mukherieel 
(|2004h : 

iLiddle et al ] (|2006bh . In particular we compare 
the luminosity distance at redshift z of each model given 
by 



d L (z) = c(l + z) 



r dz ' 

Jo H(z') 



(11) 



with the S N-Ia lumino s ity d istances from the latest 
catalogue of iRiess et all (|200/t l. This includes 182 SN- 
Ia, "flux- av e raged " with a Az = 0.05 binning as in 
I Wang et al.l (|2005l ) to reduce possible systematic effects 
from weak lensing. We also consider the new 57 su- 
per novae_mmin^_^om the E SSENCE Supernova Suvey 
of IWood-Vasev et all (2007) in combination with the 
38 nearby s upern ovae (with < 0.023 < z < 0.15) of 
IRiess et all (|2007h . We also consider the CMB shift pa- 
rameter R measured by the three-year W MAP experi- 
ment, R = 1.70 ± 0.03 lSpergel et all (|2006h . in combina- 
tion with the BAO measurement of the distance param- 
eter at r edshift z 0.35 dyjz = 0.35) = 1.300 ± 0.088 
Gpc fsee lEisenstein et"aD (|20Q5h ). The shift parameter R 
is defined as 

1 r z CMB Ay 1 

R = nl L hop)' (12) 

where zcmb = 1089 is the redshift of recombination. For 
the BAO measurement, the distance parameter is: 



'{zbao) 



2 (zbao) 



CZBAO 



(13) 



H(zbao) _ 

where r(z) is the comoving distance at redshift z and 
zbao = 0.35. We decide not to use data coming from 
weak gravitational lensing and galaxy clustering because 
they have the largest systematics and we prefer to be 
conservative in our analysis. 

2.3 A new algorithm for computing Bayesian 
Evidence 

Given a set of cosmological data, we evaluate the Bayesian 
Evidence by integrating the likelihood distribution with 
a method based on a m odified v ersion of the VEGAS al- 
gorithm. Introduced bv lLepagel {l978), VEGAS is widely 
used for multidimensional problems which occur in el- 
ementary particle physics. VEGAS is an importance- 
sampling algorithm, where regions where the integrand 
has large absolute value are sampled with a higher den- 
sity of points than regions where it is low. The key ele- 
ment of the VEGAS algorithm is that samples are drawn 
from a probability distribution which is separable in the 
coordinates. This reduces complexity in two ways. In a 
d- dimensional problem, the probability distribution p is 
specified by d one- dimensional distributions, rather than 
one d- dimensional distribution: 



p(e)(xg 1 (9 1 )g 2 (02)...9d(0d). 



(14) 



Secondly, the generation of the samples is simplified by 
successively drawing from each of the d probability dis- 
trib utions. The o ptimal weight function can be shown to 
be (|Lepagdll978h 



/ d d ~ 



PiP) 



(15) 



where f{0) is the integrand to be sampled. An initial 
(e.g. random) sampling estimates / roughly, from which 
an initial set gi can be constructed. Subsequent samples 
drawn from the q i can then be used to refine the g%. See 
IPress et al.l (1986) for further details. A problem of VE- 
GAS is that it may not do well when the integrand is 
concentrated in one-dimensional (or higher) curved tra- 
jectories (or hypersurfaces) , unless these happen to be 
oriented close to the coordinate directions. 

To solve this problem we have generalised the algo- 
rithm, and calculate, for a first and preliminary sampling, 
the covariance matrix of the Likelihood function. This we 
diagonalise and use the eigenvectors to define new pa- 
rameters. In terms of these new parameters, the likelihood 
should be closer to separable. Note that this modified VE- 
GAS algorithm should be particularly efficient if the like- 
lihood is single-peaked. For multimodal likelihoods with 
more significant peaks, it will become increasingly less ef- 
ficient, depending on the number, shape, relative height 
and orientation of the peaks. 

Given a cosmological dataset and a theoretical frame- 
work, the Likelihood function will be clearly a function of 



L(9) oc exp 



x 2 (0) 



and the algorithm calculates 
1 = j dOL{0). 



(16) 



(17) 



We always consider flat priors in our analysis; this implies 
that the priors will be constant so the Evidence will be 
given by: 



P- J, 



(18) 



where P is simply the product of the various priors for 
the parameters, P — nf=i Pi- 

In the following we describe the principal steps of our 
algorithm: 

• We do a first sampling of the likelihood function with 
Ncov sampling points and we calculate the covariance 
matrix C — {(0 — 0){0 — 6) T ), which is symmetric and 
positive definite 

• A square matrix which is symmetric and positive def- 
inite can be written as the product of a lower triangular 
and an upper triangular matrices (Cholesky decomposi- 
tion) : 

C = Q-Q T . (19) 

We use the CHOLDC routine fse dPress et all (jl986l ). 
par. 2. 9) to calculate the matrix Q. In general we can write 
the likelihood function as: 



L0) c< exp f-lfftp- 1 ] 



(20) 



where C is the covariance matrix; if y — Q 1 we now 
have: 
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NUMBER OF ITERATIONS 



NUMBER OF ITERATIONS M=5 




3000 4000 5000 6000 7000 8000 

NUMBER OF SAMPLING POINTS 



9000 1 0000 



Figure 1. Standard deviation in function of the number of iter- 
ations m used (with N = 5000 for each iteration) in the calcula- 
tion of the Likelihood function with 3 parameters (f2 m , w. Ho). 
The value of the integral is: I = 0.04915 

L(y) oc exp ^-Ijf^ ( 21 ) 

• We now choose to sample y rather than so our 
sampling should be very efficient. In fact, most sampling 
points will be generated in the subspace of the parameter 
space where the likelihood function is not zero. 

It is clear that, changing our variables, the integral of 
the Likelihood will be given by: 

/ L(0)d0 = f L(y)|det(Q- 1 )|dy. (22) 

J D J D' 

We perform m statistically independent evaluations of the 
new function using N sampling points for each iteration. 
The iterations are independent but they do assist each 
other because the algorithm uses each one to improve the 
sampling grid for the next one. The results of m iterations 
are combined into a single best answer and its estimated 
error, by standard inverse variance weighting. We also 
compute x 2 to check that the best-fitting solutions are 
acceptable statistical fits. 

Our results are clearly dependent on two parame- 
ters, the number m of iterations and the number N of 
sampling points for each iteration. The more iterations or 
sampling points are used, the more accuracy is reached. 
In Figures [T]|2] we show how much the standard devia- 
tion depends on m and N. In Figure [3] we can see the 
density of sampling points in the Q m — Ho plane; the den- 
sity will be greater in the subspace where the Likelihood 
function lives and we can also see the principal directions 
of the function. In general, implementing the routine for 
the calculation of the covariance matrix, the relative er- 
ror ^j- in the calculation of the integral is lowered by a 
factor 10, for the same number of function evaluations. 
In our analysis we always use m — 5 iterations and from 
N — 5000 to 30000 sampling points (depending on the 
dimension of the parameter space and on the size of the 
prior space), except for MODEL V for which N = 100000. 
We reach an uncertainty of ~ 10 -3 in A \n(E). We also 
used Ncov = 20000 in the first iteration for calculating 
the covariance matrix: this large number is justified by 



Figure 2. The same of Figure [T] but the standard deviation is 
plotted in function of the number N of sampling points (with 
m = 5) 



DENSITY OF SAMPLING POINTS 




Figure 3. Density of sampling points in the 2-dimensional 
space Qm — Ho 

the fact that it's important to achieve a good estimate 
of the covariance matrix, to have good "principal direc- 
tions" for the next samplings, especially when we handle 
with several parameters. 



3 RESULTS 

Let us first analyze the full iRiess etall (|2007h data in 
combination with CMB and BAO. The main results of 
this analysis are reported in Table 1 and in Figures 4-8. 
The best-fitting parameter values are the means obtained 
from the full posterior probability distribution, rather 
than the maximum likelihood values. The standard de- 
viations are similarly obtained from integration over the 
posterior. 



As we can see a cosmological constant is preferred 
by the data and is always compatible with it, indepen- 
dently of the model considered. The constraints we ob- 
tain on the equation of state parameter, assumed as a 
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Table 1. The parameter constraints together with the mean value of A \n( E) and the minimu m chisquared 
for the six models considered using the 182 supernovae by Riess et al. 2007 Rie ss et al ] (l2007h binned with 
Az = 0.05 up to z = 1.7 (about 34 bins). The (unnormalized) ln(E) for the ACDM model is \n(E) = 
— 16.2466 ± 0.0002. The uncertainty in the value of A \n(E) is calculated making use of the usual formula 

for the propagation of the uncertainty of a variable x = x(y) — > a x = ^J^^L-^ ^^^-^ <r%.. In our case 
A H E) = ln(f ) CTA ln(B) = ^) 2 + (^) 2 . 



Constraints 


AlnE 


v 2 


Model 


Q m = 0.28 ±0.03 
H = 64.5 ± 0.09 


0.0 


24.39 


I 


n m = 0.27 ±0.03 
H = 63.4 ± 1.1 
w < -0.84 at lcr 
w < -0.73 at 2(7 


-0.222 ±0.005 


22.43 


II 


Q m = 0.27 ±0.03 
Ho = 63.4 ± 1.1 
w = -0.86 ±0.1 


-1.027 ±0.002 


22.43 


III 



Qm = 0.28 ±0.04 -1.118 ±0.015 21.47 IV 

Ho = 63.8 ± 1.4 
wo = -1.03 ±0.25 
w a = 0.76^0^ 



Qm = 0.27 ±0.03 -1.059 ±0.008 21.38 V 

Ho = 63.5 ± 1.1 

WQ = -0.85 ±0.12 

wi = -0.81 ±0.21 

a s unconstrained 

q unconstrained 



n m = 0.30 ±0.05 
Ho = 63. 5+ 
wo = 

wi = 0.78^ 



-1.834 ±0.006 21.52 



VI 



-1.2 

-1.08+ - 24 



0.30 
+0.83 
0.57 



72 
70 

68 

~ 66 
o 

_S 64 
J 62 
60 

56 
0. 




Model I 
Riess07+WMAP+BAO 



10 0.15 0,20 0.25 0.30 0.35 0.40 0.45 0.50 




0,10 0.15 0.20 0.25 0.30 0.35 0,40 0.45 0.50 



Figure 4. Likelihood contours at 68% and 95% (two- 
parameter) for MODEL I 



Figure 5. Likelihood contours at 68% and 95% for MODEL 
II 



constant, are w < —0.84 in the case of w ^ —1 and same orde r of magnitude of th e previous constraints re- 
w = —0.86 ±0.1 at 68% c.l. when models with w < — 1 ported by Liddle et al. (2006b) but where the new SN-Ia 
are included. Those constraints are compatible and of the dataset o f lRiess et al. (2007) was not considered. However 
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Model III 

Riess07+WMAP+BAO 



-2,0 1 ' ' 1 1 1 1 1 1 1 1 ' 1 1 ' 

0,10 0,15 0.20 0.25 0.30 0.35 0,40 0,45 0.50 



Figure 6. Likelihood contours at 68% and 95% for MODEL 
III 




-2,0 -1,8 -1,6 -1,4 -1,2 -1,0 -0,8 -0,6 -0.4 -0,2 0,0 



Figure 7. Likelihood contours at 68% and 95% for MODEL 
IV 




Model VI 

Riess07+WMAP+BAO 



_K8 -1,6 -1,4 -1,2 -1,0 -0.8 -0,6 -0,4 -0.2 



Figure 8. Likelihood contours at 68% and 95% for MODEL 
VI 



the evidence for the w < — 1 case is worse by A \n(E) ~ 1, 
i.e. there is no indication from the data that we should 
extend the parameter space to these phantom dark energy 
models. 

The same happens when we consider models with an 
equation of state parameter which varies with redshift. 

For the models I-IV, which have the best Evidences 
and the same parametrization (with wo = — 1 and w a = 
for MODEL I and w a = for MODEL II-III) we have 
also considered results comi ng from th e Bayes ian model 
averaging; as explained in Liddl e~et al.l (|2QQ6bh . because 
of our ignorance about the true cosmological model, we 
may think that the probability distribution of the param- 
eters is a superposition of its distributions in different 
models, weighted by the relative model probability, as in 
quantum mechanics, where the state of a physical system 
is a superposition of its possibilities until a measurement 
determines the collapse in a single eigenstate. 

If we convert the A \n(E) into posterior probabilities, 
assuming equal prior probabilities, we have 40.2%, 32.2%, 
14.4%, 13.2% for models I-II-III-IV. 

The constraints on the cosmological parameters from 
the Bayesian model averaging of models I-II-III-IV are: 
On = 0.27 ± 0.03, Ho = 63.2+^, w = -1.0 ± 0.1 at la, 
w a = 0.01q q3 at 2cr. The confidence limits in w a are exa- 
cly zero at la; this is because the probability distribution 
for this parameter is a delta function for models I-II-III 
and it is superimp osed to the extend ed tails of model IV 
(as explained in lLiddle et al.l (|2006bh ). 

It is also interesting to consider the effects on the 
cosmological parameters if we remove from the dataset 
the supernovae with z > 1. In Table 2 we report our 
principal results. 

As we can see, there are no significant differencies 
on the mean values however the error bars are generally 
reduced by a ^ 30%. Models with varying- with redshift 
equation of state have a slightly better evidence but with 
a cosmological constant is still favoured. 

It's also useful to check if our results are the same 
when we consider a different dataset. To this extent, we 
use the 57 sup ernovae coming from the ESSENCE Su- 
pernova Suvey IWood -Vasev et al. in combination 
with the 38 near by supernovae (wit h 0.023 < z < 0.15) of 
Riess et al. 2007 (jRiess et al.l (|2Q07h ). We do not consider 
models with evolving dark energy, because in this analysis 
we limit our redshift range to z < 0.670. 

The results are reported in Table 3. The results are 
fully compatible with those from the previous analysis 
and, again, they provide a substantial evidence for a con- 
stant w, with a cosmological constant being preferred. 



4 CONCLUSIONS 

In this paper we have carried out a Bayesian model 
selection analysis of several dark energy models using 
the n ew data of h i gh re dshift supernovae of Riess et al. 
2007 (Rie ss et al.l (120071 ) ] and from the ESSENCE sur- 
ll (L 



vey (| Wood- Vasev et al.l (|2007l )). together with Baryonic 
Acoustic Oscillations and Cosmic Microwave Background 
Anisotropics. To this extent, we have developed a new al- 
gorithm to calculate the Bayesian Evidence which is fast 
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Table 2. The parameter constraints together with the mean value of A \n(E) and the minimum chisquared 
for the six models considered when z < 1, using the 166 supernovae with z < 1 of Riess et al. 2007 binned 
with Az = 0.05 (about 20 bins). 



Constraints 


A InE 




Model 


Qm = 0.27 ±0.04 


0.0 


16.43 


I 


li q — 04. o zt u.uy 








Qm = 0.27 ±0.04 


-0.223 ±0.004 


14.54 


II 


H = 63.4 ± 1.3 








w < -0.82 at la 








w < —0.73 at 2a 








Qm = 0.27 ±0.04 


-1.022 ±0.004 


14.54 


III 


H = 63.4 ± 1.2 








w = -0.86 ±0.1 








Qm = 0.29 ± 0.04 


-1.090 ±0.011 


13.33 


IV 


H = 63.6 ± 1.2 
wo = -0.98±g;|| 
w a = 0.72^-- 








Qm = 0.28 ±0.04 


-1.020 ±0.008 


13.00 


V 


H = 63.4 ± 1.0 








wo = -0.88 ±0.12 
W! = -0.63l°-j 7 
a s = (unconstrained) 
q = (unconstrained) 








Q m = 0.30 ±0.05 


-1.691 ±0.01 


12.95 


VI 



Ho = 64.3 ± 1.2 
w = -1.17 ±0.37 
w\ = (unconstrained) 



Table 3. Paramete r constra i nts a nd Evidence for MODEL I-II-III, using the 95 supernovae of 
ESSENCE+nearbv [Riess et al.l (|2007l ) binned with A z = 0.05 up to z = 0.65 (15 bins) 



Constraints 


A InE 


X-Min 


Model 


Q m = 0.25 ±0.04 
Hq = 64.7 ±0.09 


0.0 


10.73 


I 


Q m = 0.27 ±0.04 
Ho = 63.6 ±1.3 
w < -0.80 at la 
w < -0.67 at 2a 


-0.258 ± 0.004 


9.28 


II 


Qm = 0.27 ±0.04 
Ho = 63.5 ±1.2 
w = -0.86 ±0.11 


-1.027 ±0.003 


9.28 


III 



(less than 1 hour for each calculation of the Evidence on 
a 2GHz CPU) and very accurate (a relative uncertainty 
^ ln(E) < 10~ 2 in 10 5 likelihood evaluations). 

We find that with current observational data the 
usual ACDM model is slightly preferred with respect to 
dark energy models with equations of state in the range 
— 1 ^ w ^ and substantially preferred to dark energy 
models with — 2 ^ w ^ or to dark energy models with 
an equation of state which evolves with time. 

However we would like also to stress that it may be 
premature to reject models only on the basis of Bayesian 
model selection; in general, the simplest models may not 



be the "true" model. In this way, until we have a the- 
oretical explanation of the accelerated expansion of the 
Universe, one should keep an open mind to all the alter- 
natives to the ACDM scenario, even if, at the moment, it 
seems the simplest description of our Universe. 
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